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Abstract 



The collimation of average multiplicity inside quark and gluon jets is investigated in per- 
(~| \ turbative QCD in the modified leading logarithmic approximation (MLLA). The role of 

higher order corrections accounting for energy conservation and the running of the cou- 
pling constant leads to smaller multiplicity collimation as compared to leading logarithmic 
approximation (LLA) results. The collimation of jets produced in heavy-ion collisions has 
' also been explored by using medium-modified splitting functions enhanced in the infrared 

■ sector. As compared to elementary collisions, the angular distribution of the jet multiplicity 
^ ! is found to broaden in QCD media at all energy scales. 
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■ 1 Introduction 



Paradoxically, the "jet quenching" phenomenon observed at RHIC precedes by 
many years the first measurements of jets in heavy-ion collisions. The spectacular 
suppression of large-p^ single-inclusive pion production in central Au-Au colli- 
sions has been widely interpreted as coming from the energy loss of hard partons 
in a dense medium (see e.g. [1]). However, these data - as well as more differen- 
tial measurements such as di-hadron and photon-hadron correlations - only hardly 
inform us on quark and gluon multiple scattering processes and medium-induced 
gluon radiation in quark-gluon plasma (for a review on the subject, we refer the 
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reader to [2]). Measuring hadronic distributions inside reconstructed jets in heavy- 
ion collisions should therefore be key in order to access the underlying dynamics 
of jet quenching. Using recent advances in jet algorithm techniques [3], prelim- 
inary measurements on inclusive jet spectra have been reported lately in Au-Au 
collisions at y/s^ = 200 GeV by the STAR collaboration [4]. On the theoretical 
side, various jet observables have been investigated such as the angular distribu- 
tion of jet average multiplicities [5,6], multiplicity distributions [7,8] or inclusive 
momentum spectra [9]. The appearance of various parton showers in heavy-ion 
collisions [10-12] should also allow for a more systematic exploration of jet ob- 
servables and their medium- modifications. 

As for any hard process in heavy-ion collisions, the precise knowledge of the ex- 
pected baseline in p-p collisions becomes crucial in order to properly quantify the 
effects of quark-gluon plasma formation on jet physics. Over the past twenty years, 
significant progresses have been achieved in order to improve predictions on vari- 
ous jet observables. In particular, the modified leading logarithmic approximation 
(MLLA) in perturbative QCD has been successfully tested from e+e^ to hadronic 
collision^ [14]. 

To this sake, we investigate in this paper the coUimation of the average multiplicity 
inside a quark and a gluon jet as a function of the jet energy scale. The coUimation 
is characterized by the cone aperture of the sub-jet 65 that contains a fraction S 
of the jet average multiplicity. Let N{E, Bq) be the multiplicity in a jet of energy 
E and opening angle 60, and iV(6; E, Go) that of the sub-jet of opening angle 6 
inside this jet, Qs is determined by solving the multiplicity coUimation equation, 

N{Qs;E,Qo) = SxN{E,Qo)- 

In the leading logarithmic approximation (LLA) [15], the solid angle 9 1/2 con- 
taining half of the average multiplicity of the jet 60 decreases with the jet hard- 
ness Q ~ EQo approximately like 0i/2(<5) ~ N^^/^(Q), such that at high en- 
ergy scales the bulk of the total multiplicity is concentrated at smaller solid angles 
around the direction of propagation of the jet. In this paper, we extend this cal- 
culation to the MLLA by incorporating all corrections of order O (^y/afj, which 
partially guarantee energy conservation and account for the running of the cou- 
pling constant in intra-jet cascades. The treatment of such a complicated task 
in perturbative QCD simplifies tremendously thanks to the angular ordering of the 
successive emission of soft gluons, which leads to simplified jet evolution equa- 
tions [15]. At the end of the cascading process, the coUinear cut-off parameter Qo 
can be taken as low as A^^^^^ (the so-called limiting spectrum approximation), and 
the local parton-hadron duality (LPHD) hypothesis [16] can then be advocated so 
as to determine hadronic distributions. 

^ For rather exclusive observables such as the transverse momentum of hadrons inside jets, 
the disagreement between data and MLLA expectations can be cured by the inclusion of 
higher-order terms, O (a^), in the calculation [13]. 
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Within the same fomialism, the coUimation of jets has also been detennined in 
heavy-ion collisions by inserting medium-modified splitting functions, in which the 
soft gluon emission is arbitrarily enhanced from a toy QCD-inspired model [9]. Al- 
though more realistic approaches treating parton energy loss in medium-modified 
fragmentation have been proposed (see e.g. [17] for a review), the use of this model 
was motivated by the fact that analytic calculations can more easily be performed. 

The outline of the paper is as follows. In Section 2, the calculation of the average 
multiplicity in a sub-jet N is carried out at MLLA, which eventually allows for 
solving numerically the above coUimation equation. In Section 3, our results are 
compared to the LLA predictions in the vacuum and medium effects on jet coUi- 
mation are discussed. 



2 Theoretical framework 



2. 1 Energy-multiplicity correlation and sub-jet average multiplicity 



Consider a jet produced in a high energy collision, initiated by a parton of flavour 
Aq and energy with an opening angle 6o which separates it from other jets. 
Inside this jet, a sub-jet initiated by a parton A and energy uE is defined by the 
opening angle © < ©o- The sub-jet multiplicity TV^^, i.e. the mean number of 
hadrons produced inside the angular range © < ©o of the jet Aq, is given by [15, 18] 

7Vi(©;£;,©o)« E /' duuDX{u,EQo,uEQ)N\{uEQ,Qo) , (1) 



'Qo/EO 

which determines the correlation between the jet axis and the average multiplicity 
of the sub-jet. In Eq. (1), D^^ is the probability to find a parton A with energy 
fraction u and virtuality uEQ > Qq inside the jet initiated by the parton Aq, and 
N\{uEQ, Qq) is the bare average multiplicity in a sub-jet of energy uE and open- 
ing angle 0. Since soft particles are less sensible to the energy balance in intra-jet 
cascades, the correlation between the energy flux and the sub-jet multiplicity dis- 
appears. As a consequence, the measured average multiplicity iV^^ factorizes in the 
form [18] 

<(leo,le) ^ ^{C),^iYe,,Ye)N^iYe), (2) 

c 

after the Taylor expansion of D^^ in Eq. (1), where we define 

EQo ^ ^ EQ 
>Go = In ^r— , Ye = In — . 

In Eq. (2), {C)ao is the colour current of partons caught by the calorimeter and 
Ng{YQ) is the bare average multiplicity of a gluon jet of energy E and opening 
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angle 6, given by the solution of the MLLA evolution equations. The colour current 
{C)ao describes the evolution of the jet between the scales Qeo = EQq and Qq = 
EQ and indicates that the registered partons lose the memory of the initial colour 
state of the parent parton Aq because of intra-jet evolution [15,18]. Introducing the 
gluon to quark jet total multiplicity, r = Ng /N^, the colour current can be written 
at MLLA accuracy as [18] 



d In 



dY, 







where for the sake of brevity, we introduced the variable[il 



1 ,jye. 



(3) 



(4) 



with A^^ = 3 is the number of colours, (3q = ^ (y^c ~ f j the first coeffi- 
cient of the QCD /3-function and Tjj = n//2 where n/ = 3 is the number of active 
flavours. The functions appearing in Eq. (3) are written in the form 



(0 K iu,o = (j = 2, , 



i=2 



(5a) 
(5b) 



which can be determined from the DGLAP evolution equations [18]. 



2.1.1 Medium-modified DGLAP and MLLA evolution 

In order to compute Eq. (3) from Eq. (5), the moments V{j,^) need first to be 
determined by solving the DGLAP equation in Mellin space. 



d 



-pUj) 

\ V,,{j) VggU) ) 



(6) 



where Vq^^ and Vq^ stand respectively for the flavour non-singlet (or valence) and 
flavour- singlet quark distributions, and Vik{j) is the Mellin transform of the leading- 
order splitting functions Pik{z). In order to account for the medium-induced gluon 
radiation in heavy-ion collisions, various attempts to determine medium-modified 



It should not be confused with the notation ^ = ln(l/x) used e.g. in [19]. 
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splitting functions have been recently performed [20-22]. In this paper, we shall 
adopt the most simple approach proposed by Borghini and Wiedemann [9] which 
allows for analytic solutions. In this model the infrared sensitive parts of the split- 
ting functions are enhanced by a factoid] A^^, 



z 



1 - z 



+ z(l-z) 



2N. 



z 



+ Z-2 



2N. 



VI- z 



2T^[z' + {l-zn 
(7a) 

-1-^), (7b) 



with the [...]+ prescription defined as dx[F[x)]^g{x) = dxF{x)[g{x) — 
g{X)\- Performing the Mellin transform of Eq. (7a,7b) gives [23] 



P,,(j) = -4iV^ 



N, + 1) + Ns-fE - 



3 



2nf ^8NSf+j + l) 



3 3 ■ jif-l)ij + 2) 
3^+3 + 2 



j(j + l)(j + 2)' 

(2Ns-l)if+j)+2 



r,,{3) = -c. 



4Ns + 1) + 4iV,7s - 4 



AT, -1 



3(3 + 1] 



(8a) 
(8b) 
(8c) 
(8d) 



where ip{j) is the digamma function. It can be easily checked that Eq. (8) reduces to 
the ordinary splitting functions given in [15,24] after setting Ng = 1. Using Eq. (8), 
the DGLAP equation Eq. (6) can be solved, such that the expressions for V^(j,^) = 
1^gj3,0 + ^^fUO, '^1(3,0 and Vl{j,0 can be taken from [15] by 

introducing the A^^ dependence from Eq. (8). Finally, these distributions and their 
derivatives can be evaluated at j = 2 in order to determine the moments (u) and 
{u \nu), cf. Eq. (5). 



The other ingredient that enters the evaluation of the colour currents Eq. (3) are the 
MLLA gluon and quark jet average multiplicities [15, 24, 25]. In the model of [9], 
they read [26] 



dY, 



e 



/C Fq-'^/^o exp 
A^. 



Ye 
(3o ' 



CI. 7o 



(9a) 
(9b) 



In [9] the parameter /^^^ is used, corresponding to A'^s — 1 here. 
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with 0"! ^ 0.28, and ri ~ 0.185 for nj = 3. The anomalous dimension, 70, deter- 
mines the rate of the multiplicity increase in a jet; it can be written 

,HEe) = 2K = (10) 

TT (Jo Ye 

The formulae in Eqs. (9a,9b) follow from the MLLA evolution equations as a con- 
sequence of angular ordering in the partonic shower. The constant /C normalizes 
the number of partons to the number of charged hadrons according to the LPHD 
hypothesis [16]; as can be checked below, this constant does not play any role in 
this context. Finally, working out the structure of Eq. (3), one has 



(11) 



+ 70 X 



which will be used in the present form in the following. As expected, the MLLA 
expressions of the colour current in the vacuum [18] are recovered when setting 
A^, = 1 inEq. (11). 



2.2 Collimation of average multiplicity inside a jet 



As stressed in the Introduction, the collimation is characterized by the angular size 
of the cone containing the fraction 5 < \ of the total multiplicity in a jet, 

N\,{Ye,,Ye,)=5xNX^{Ye,). (12) 

Using Eq. (2), the collimation equation (12) becomes 

(C),^(yeo, VeJiV^'(leJ = 5 x NX,{Ye,), (13) 

c 

which can be solved numerically by using the analytic expressions for the colour 
currents and the total jet average multiplicity. 

The collimation of average multiplicity in a jet has been determined in [15] at 
At large opening angles ©o ~ 1, the LLA collimation was found to scale 

like 

[iV^(E/A,,j]^^*''^% (14) 

indicating that the bulk of the jet multiplicity concentrates closer to the jet axis as 
the energy scale increases. For the sake of more accurate predictions, this result is 
extended in the present paper to MLLA accuracy, i.e. keeping track of all terms of 



Note that the collimation of the energy has also been studied at LLA in [15]. 
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order O (^v^olj in (13). Moreover, by inserting an explicit A^^ -dependence in the 
splitting functions, we provide results accounting for enhanced soft gluon radiation 
in dense QCD media (Ng > 1). 



3 Phenomenology 

3.1 Collimation from LLA to MLLA 

In this Section, the multiplicity collimation of jets produced in the vacuum is com- 
puted at MLLA and compared to the LLA prediction, recovered by setting cti, ri 
and 7o to zero in Eq. (13). In Fig. 1, the normalized angular aperture, ©^/©o, is 
plotted as a function of the jet energy scale Q = EQq, in units of Aq^^^, for gluon 
(left) and quark (right) jets at LLA (solid lines) and MLLA (dashed) accuracy. For 
the illustration, ©^ has been determined using 5 = 0.25 and 0.5. 




Figure 1. Collimation of average multiplicity inside a gluon jet (left) and a quark jet (right) 
in the vacuum at LLA (solid) and MLLA (dashed) accuracy. 

Although the trends are similar, ©^/©o turns out to be larger at MLLA than at LLA, 
indicating that corrections following from energy conservation and the running of 
the coupling constant increase the size of the angular distribution, i.e. making jets 
less coUimated at MLLA. As expected, MLLA corrections are larger at low jet 
energy scales and closer to LLA at asymptotic energies as the anomalous dimension 
vanishes. Also, MLLA corrections prove slightly larger for a gluon than for quark 
jet, yet the difference is rather small. More remarkably, the difference between LLA 
and MLLA predictions is somewhat larger for smaller values of the parameter 6. 
At low scales, Q/Aq^d = ^0^' ®c is increased by 50% from LLA to MLLA for 
6 = 0.25, but only by 15% when 6 is set to 0.5. 
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These predictions could be tested at the Tevatron and the LHC, by measuring jets 
of transverse momenta say ~ 100 GeV and angular aperture 60 ~ 0.5, corre- 
sponding to energy scales Q/Aq^^ ~ 200. Lower scales, for which MLLA predic- 
tions are largest, could even be reached at RHIC energy. 

2.2 Medium ejfects 

Let us now discuss how the multiplicity coUimation of jets produced in heavy- 
ion collisions could be distorted by quark-gluon plasma formation. Preliminary 
predictions have been made by enhancing the infrared gluon emission by 60% and 
80% (respectively, Ng = 1.6 and Ng = 1.8) and compared to the MLLA vacuum 
expectations {Ns = 1). The values of A^^ are somehow arbitrary; these numbers 
were chosen in [9] in order to reproduce the magnitude of single-inclusive pion 
suppression data in Au-Au collisions at RHIC. 




Figure 2. CoUimation of average multiplicity for 6 = 0.5 inside a gluon jet (left) and a 
quark jet (right) produced in the vacuum (Ng = 1) and in the medium {Ng = 1.6 and 
Ns = 1.8). 

The predictions are displayed for gluon (jets) and quark (right) jets in Fig. 2 by 
setting 5 = 0.5. As expected from the role of soft gluon radiation in jet cascades, 
the angular size Qs broadens as compared to the vacuum case. The typical values 
for the angular aperture ©^/©o is increased by roughly factor 0.15-0.2 as compared 
to jets produced in the vacuum, depending on the values of Q/A^^^ and Ng. This 
observable could be measured in heavy-ion collisions at RHIC and at the LHC 
provided the background from the underlying event can be safely removed from 
the jet signal, using the techniques discussed in [3]. On top of an increased of 
the total multiplicity, see Eq. (9a), jets produced in heavy-ion collisions appear 
to be less coUimated than in the vacuum, according to the model [9] used here. 
This observation seems qualitatively consistent with the /Cj^ -broadening of the jet 
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multiplicity distributions studied in [6], yet this jet shape is directed related to the 
coUimation of energy rather than to the coUimation of multiplicity considered here. 

The calculations of the medium-modified coUimation of multiplicity discussed in 
this Section are not meant to be quantitative, given the rather primitive model used 
to describe the process of parton multiple scattering and induced-gluon radiation. 
Nevertheless, the qualitative prediction of a reduced coUimation of jets produced 
in heavy-ion collisions is expected to hold in more realistic calculations. The pro- 
cedure presented in this paper in order to compute the coUimation of multiplicity 
at MLLA from given splitting functions could for instance be used in the future to 
provide more accurate predictions in heavy-ion collisions. 



4 Summary 

In this paper we have considered the coUimation of multiplicity distributions inside 
jets produced at high energy colliders. The calculation performed at LLA in [15] 
has been extended to MLLA accuracy by taking full care of all terms of order 
O (^v^) evolution. The sub-jet containing a fraction 6 of the total jet 

multiplicity has been found to widen at MLLA as compared to previous results at 
LLA. The difference between LLA and MLLA expectations are largest at small jet 
energy scales and could be tested at present (RHIC, Tevatron) and future (LHC) 
hadronic colliders. 

The same analysis was performed by using the medium-modified splitting funtions 
introduced in [9], which enhance the soft gluon radiation of a highly virtual parton 
traveling through quark-gluon plasma. The calculation has also been performed at 
MLLA and compared to the vacuum expectations. Results indicate that jets pro- 
duced in heavy-ion collisions are expected to be less coUimated than in the vac- 
uum, even though quantitative predictions using more realistic medium-modified 
splitting functions are yet to be performed. 
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